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ABSTRACT 

We revisit the formation and evolution of the first galaxies using new hydrodynamic cosmological simulations 
with the ART code. Our simulations feature a recently developed model for H2 formation and dissociation, and 
a star formation recipe that is based on molecular rather than atomic gas. Here, we develop and implement a 
new recipe for the formation of metal-free Population III stars. We find the epoch during which Pop III stars 
dominated the energy and metal budget of the first galaxies to be short-lived. Galaxies which host Pop III stars 
do not retain dynamical signatures of their thermal and radiative feedback for more than 10'^ yr after the lives of 
the stars end in pair-instability supernovae, even when we consider the maximum reasonable efficiency of the 
feedback. Though metals ejected by the supernovae can travel well beyond the virial radius of the host galaxy, 
they will typically begin to fall back quickly, and do not enrich a large fraction of the intergalactic medium. 
Galaxies more massive than 3 x 10^ Mq re-accrete most of their baryons and transition to metal-enriched Pop 
II star formation. 

Keywords: galaxies: formation — galaxies: evolution — methods: numerical — stars: formation — cosmol- 
ogy: theory 



1. INTRODUCTION 

A natural consequence of Big Bang Nucleonsynthesis is 
that the first stars in the universe formed from gas that was 
completely devoid of elements heavier than lithium. Since 
such conditions do not exist in any known star-forming re- 
gions today, star formation in the primordial regime has thus 
far only been explored by theory and simulations. From first 
principles, it can be deduced that gas which is free of met- 
als would not be able to cool efficiently, and would there- 
fore inherently have a higher Jeans mass than stars forming in 
metal-enriched gas. iBromm et alj (119991) showed that metal- 
free gas settles into disks, then fragments into clumps with 
Mj « IO^Mq, which undergo runaway collapse to densities 
of riH > 10^ cm"^. In a follow-up study, it was shown that 
this process was robust to initial conditio ns in Smoothed Par - 
ticle Hydrodynamics (SPH) simulations (iBromm et al.ll2002l) . 
Using indep endent Ada ptive Mesh Refinement (AMR) grid 
techniques, lAbel et all ([2000) confirmed that such dense 
clumps and cold pockets can indeed form in primordial gas. 
lAbe l et all (2002) presented even more realistic simulations, 
and concluded that single massive stars would form via H2 
cooling at the center of these clumps, and the radiative feed- 
back would halt accretion onto the star and prevent further star 
formation in the parent halo. However, recent studies with 
both SPH (Clark et al. 20 U) and grid techniques (Greif et al. 
I2OI ItlHosokawa et al.ll201ll) have suggested that angular mo- 
mentum imparted on the gas during collapse could still lead 
to fragmentation, causing the cores and the resulting stars 
to be substantially smaller. Ultimately, the next generation 
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of super-zoomed simulations will need to follow the proto- 
stellar systems for longer periods of time, with more detailed 
treatments of radiative transfer and magnetohydrodynamics, 
to conclude decisively on the true nature of Pop III stars 
fGreifet al. 2012). 

Though we do not yet have any direct observational evi- 
dence for Pop III stars, we know that they must have existed in 
some form, as gas in the universe inevitably transitioned from 
having primordial composition to being enriched with enough 
metals to allow present-day star formation to commence. The 
nature of this transition is of key importance for understand- 
ing the dawn of galaxy formation: if Pop III stars did indeed 
form with a top-heavy initial mass function (IMF), a signifi- 
cant fraction of them may have been prone to end their lives 
as pair-instability supernovae (PISNe). Such supernovae gen- 
erate up to ten times as much thermal energy as Type II SNe 
(jHeger & Wooslev 2002), quickly heating the gas in their host 
halos. In addition, metal-free stars are able to produce enor- 
mous quantities of ionizing photons: a metal-free star will 
always have a higher sur face temperature than an enriched 
count erpart of equal mass ( lSchaererf2002l : lTumUnson & Shulll 
I2OOOI) . If the IMF is indeed top-heavy, the effect is even more 
drastic as many stars would have their emission spectra peak 
in the UV regime. Both the supernovae and the ionizing pho- 
tons serve to heat and disperse the gas surrounding the region 
where the star forms, delaying the onset of steady, continu- 
ous star formation. On the other hand, PISNe release a large 
amount of metals, rapidly enriching previously primordial gas 
(Wise & Abel 2008; Greif et al. 2010; Wise et al. 2012), and 
potentially lead ing to quick termination of the Pop III epoch 
dYoshida et al.l 12004) . The balance of these effects works to 
determine star formation rates in galaxies and their subse- 
quent evolution. 

While Pop III stars are no longer thought to be a major 
driver of the global reionization of the intergalactic medium 
dCiardiet al.1 120001: fecotti et al. 2002; Ricotti & Ostrike^ 
2004; Mesinger et al. 2009), constraints on the cosmological 
electron scattering optical depth from WMAP suggest that ha- 
los less massive than 10** M© may have contributed to the 
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photon budget at the beginning of reionization. In order for 
this scenario to work, low mas s halos must permit the escape 
of ion izing photons effectively (lAlvarez et alj2012tlAhn etaTI 
l2012h . Pop III stellar feedback has been explored as a mech- 
anism to create windows of time during which such high es- 
cape fractions are made possible (Wise & Cen 2009 , but see 
iGnedin et al.ll200l iRicotti et al.ll2008lK 

We revisit these important conclusions with novel high- 
resolution cosmological simulations that feature separate star 
formation criteria for Pop III and Pop II stars. While we can- 
not resolve all stages of their formation, we do resolve the 
clumps of gas on a ~ 1 pc scale which inevitably collapse into 
Pop III stars. In addition to the commonly accepted parame- 
ters for Pop III formation and feedback, we also explore mod- 
els in which the radiative and SNe feedback are taken at ex- 
treme values, as well as models in which the IMF of Pop III 
stars is taken to be identical to Pop II counterparts. 

Our analysis in this paper focuses on the dynamical effects 
that Pop III stars can impart onto their host galaxies. We quan- 
tify the ability of Pop III stars to suppress star formation, ex- 
pel gas, and enrich the medium both within and outside of 
the galaxies in which they form. Using a wide suite of sim- 
ulations, we show that mass resolution and mesh refinement 
criteria affect the derived importance of Pop III stars. In a 
forthcoming second paper (Muratov et al. 2012b, Paper II), 
we will explore the nature of the transition between Pop III 
and Pop II star formation, and assess the relative importance 
of feedback effects from the two stellar populations. 

2. SIMULATIONS 

We perform cosmological simulations with the Eule- 
rian, gasdy namics-HN-body adaptive refinement tree (ART) 
code (Kravtsov et al.l[T997t lKravtsovl[T999l l200l iRudd et alJ 
l2008h . The latest version of the code incorporates a new phe- 
nomenological prescription for molecular hydrogen forma- 
tion on dust grains and self-shie lding, as well as shielding by 
d ust in troduced in Gnedin et al. (2009) and developed further 
in lGnedin & Kravtsov (201 1) . Having such a detailed account 
of molecular gas, as well as excellent spatial resolution at high 
redshift makes it practical to consider a star formation recipe 
that is also based on molecular gas. This formulation has pre- 
viously been shown to much better reproduce the Kennicutt- 
Schmidt relation for high-redshift, low-mass, and low metal- 
licity galaxies, and it now enables more realistic simulations 
of the early universe jTassis et"an 120081 : iGnedin & Kravtsovl 
120101 120T1I) . Radiative transfer, including Lyman- Werner H2 
dissociating feedback, is computed using the OTVET approx- 
imation (Gnedin & Abel 2001). 

For our baseline runs, we use a 1 /;"' Mpc comoving box 
and the WMAP-7 cosmology (f7,„ = 0.28, f^A = 0.72, h = 0.7, 
CT8 = 0.817, ni, = 0.046, = 0.234). Additional runs were 
performed with 0.5 and 0.25 Mpc comoving boxes to 
explore the effects of mass resolution (hereby referred to as 
H Mpc and Q Mpc runs, respectively). The details for each 
run performed in our simulation suite are presented in Sec- 
tion 12.31 Numbers quoted in this paragraph, as well as Sec- 
tions |2T| and |22] refer to our baseline 1 /i"' Mpc runs. These 
runs start with a 256^ root grid, which sets the DM particle 
mass to moM = 5.53 x lO-^M© and the base comoving reso- 
lution of 5.56 kpc. We employ Lagrangian refinement crite- 
ria, refining cells when the DM mass approximately doubles 
compared to the initial mass in the cell, specifically, when it 
exceeds 2 x mom x x 0.8, or the gas density surpasses 



an approximately equivalent value modulated by the cosmic 
baryon fraction 0.3 x m^^vf x x 0.8. In both refinement 
conditions, the extra factor of 0.8 is the split tolerance. We 
use a maximum of 8 additional levels of refinement, giving us 
a final resolution of 10^/i~'/256/2^ « 22 comoving pc. Since 
we are studying high-redshift galaxies, it is important to note 
that this translates to about 2 physical pc at the endpoint of 
our simulations, z = 9, and 1 physical pc at z = 20 when the 
first stars begin to form. This spatial resolution is sufficient to 
capture the detailed multi-phased structure of the interstellar 
medium (ISM) (e.g. I Ceverino & Klvpinll2009l) . 

In order to simulate several representative regions of the 
universe, we employ the "DC mode" formulation presented in 
Sirko (2005) and Gnedin et al. (2011). Running simulations 
with different DC modes allows us to sample cosmologically 
over- and under-dense regions without actually changing the 
total mass within each box. A single parameter Aoc, which is 
constant at all times for a given simulation box, represents the 
fundamental scale of density fluctuations present in the box. 
At sufficiently early times, when perturbations on the funda- 
mental scale of the box are in the regime of linear growth, 
Adc is related to the overdensity. Socio) = D+{a) Aqc, where 
D+ is the linear growth factor. The expansion rate of the in- 
dividual box relates to the expansion rate of the universe by 
the fo Uowing relation, which is Equation 3 of IGnedin et alj 
(l20TTh : 

^box — ~ ~ ~ ' TTTTT'' 

[l + Aoc£'+(fl™)]'/' 

where abnx is the local scale factor of the simulation box, 
while fl„„, is the global expansion factor of the universe. For 
this study, we have used three different setups with Aoc = 
-2.57,-3.35, and 4.04, labefled 'Box UnderDense-', 'Box 
UnderDense-n' and 'Box OverDense', respectively. At the 
endpoint of our simulations, z = 9, these values translate to 
overdensities of 5dc = -0.257,-0.335, and 0.404, respec- 
tively. Boxes UnderDense- and UnderDense-n have nega- 
tive DC modes, implying they are underdense regions of the 
universe. However, while Box UnderDense- is representa- 
tive of a void, and hosts only low-mass galaxies which col- 
lapse relatively late. Box UnderDense-n hosts the first star- 
forming galaxy among all three simulation boxes. This galaxy 
is also more massive than any of those in Box OverDense until 
z« 13. 

Box OverDense hosts several massive star-forming galaxies 
which statistically dominate the sample of simulated galaxies. 
The H Mpc and Q Mpc boxes used Aoc = 5.04 and 6.11, 
respectively. These runs are primarily designed to explore the 
earliest possible epoch of Pop III star formation, tracing only 
the most overdense regions with even higher mass resolution. 
None of our simulations continue past z = 9, as the boxes are 
too small to capture the nonlinear growth of large-scale modes 
at later epochs. 

We construct catalogs of halo properties from simulation 
output s using the profiling routine described in IZemp et aP 
(I2OI2I) . We take the virial radius, 7?,.,>, as the distance from the 
center of a halo which encloses a region that has an overden- 
sity of 180 with respect to the critical density of the universe. 

2. 1 . Population II star formation 

Foll owing IGnedin et all (l2009h and IGnedin & Kravtsovl 
( I2OI ll) . we set the threshold for Pop II star formation in a 
gas cell when the fraction of molecular hydrogen exceeds 
the threshold = 2«//2/nH = 0.1. Tests described in the 
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above studies showed that the exact value of this threshold 
was not important for overall star formation rates, but mainly 
regulated the number and mass of stellar particles produced. 
The simulations performed in that study employed a non-zero 
floor (minimum value) of the dust-to-gas ratio in cells, which 
was meant to account for unresolved pre-enrichment. Since 
our current simulations spatially resolve regions where the 
first stars are expected to form, it is unnecessary and inappro- 
priate to use such a floor. This means that in our runs, prior 
to the metal feedback from the first generation of stars, only 
primordial chemistry is used in H2 formation reactions. We 
find that such primordial reactions with our resolution do not 
yield molecular fractions ///, > 0.01 on relevant timescales. 
Therefore, in order to form the first (Pop III) generation of 
stars, a separate prescription is required and is outlined in the 
next section. 

In cells where the molecular fraction exceeds the ///^ thresh- 
old. Pop II stellar particles are formed with a statistical star 
formation delay, dtsF = 10^ yr, implemented by drawing a 
random number, P, between and 1 , and forming stars only 

if F > exp {^-^^ , where dt is the length of the timestep at 

the cell's refinem ent level. Each partic le represents a stellar 
population with a'Mi ller & Scald d 19791) IMF from 0. 1 M© to 
100 Mq. The mass of a stellar particle is determined by the 
following the relation: 

P* = —PH,. (2) 

TSF 

The star formation efficiency pe r free-fall time is set to e // = 
0.01, based on recent results of iKrumholz et al.l (120121) . We 
use a constant star formation timescale Tgf = 8.4 x 10^ yr cor- 
responding to the free-fall time at hydro gen number density 
riH = 50 cm"^. This approach differs from lGnedin & Kravtsovl 
(llOll), where the timescale was computed by using the phys- 
ical density of molecular clouds. However, we find that in our 
simulations the density-dependent timescale instills a strong 
resolution dependence on the star formation rate. Pop II stel- 
lar particles are treated as statistical ensembles of stars for 
which the appropriate metal yield and fraction of stars to go 
supernovae is computed by integration of the IMF. The num- 
ber of SN II explosions is 75 per 10'* M© formed, while the 
amount of SN II metals generated by a stellar particle is 1.1% 
of its initial mass. Each supernova releases 2 x 10^' erg ther- 
mal energy which is deposited over the course of 10^ yr. 

2.2. Population III star formation 

We model the formation of Pop III stars b ased on criteria 
derived from simulations of lAbel et al.l (120021) . These authors 
showed that once the core density of a proto-cloud reached 
1000 cm"-', further collapse to a massive stellar object was 
imminent. Analyzing their results, we found that for gas at 
any given density nn past this threshold, the time of collapse 
to a stellar core is approximately six times the free-fall time 
for that density, 6fyy(n//). This collapse time is 9 Myr for 

—1/2 

nn = 1000 cm and scales as . For our fiducial runs, 
we use nfjjnin = 10000 cm"-' as the threshold and dtsF = 2.% 
Myr as the statistical star formation delay, simulating the col- 
lapse time. This value is lower than the one used for Pop II 
stars. This density threshold value was chosen to ensure Pop 
III stars would form primarily when cells have been maxi- 
mally refined, but is low enough such that the collapsing gas 
clouds are still fully resolved in our simulations. Further dis- 



cussion is presented in Section l2.3.1l 

We also set a threshold for the minimum fraction of molec- 
ular hydrogen at 10"^ to reflect that primordial gas clouds 
must cool primaril y via ro-vibrational trans i tions of H2 to 
form the first stars (Couc hman & R ees"1986l: [Tegmark et alj 
[T997I) . The precise value of this threshold is rather arbitrary, 
as we do not attempt to model the actual chemistry of stellar 
core formation. We have chosen this value because it is lower 
than, but close to the typical value for the molecular hydrogen 
fraction in cold, dense primordial gas around z = 20, which we 
have found empirically to be 2 x 10"^ (see Section |2.3.2| and 
Figure |3ll. A minimum threshold for molecular fraction en- 
sures that the //2-dissociating Lyman Werner radiation from 
recently-formed Pop III stars will realistically suppress fur- 
ther Pop III star formation in the region. 

The IMF of Pop III stars is currently a hotly debated and ac- 
tive area of research. It is still unclear whether the high Jeans 
mass of primordial gas results in a top-heavy IMF as pre- 
dicted by early sfiid ies (lAbel et al.ll2"00a iBromm et al.ll 19991; 
lYoshida et al.ll2003h . or if the angular momentum and radia- 
tive effects during infall can fr agment the cloud and gener- 
ate relatively low-mass cores (Greif et al. 201 It IStacv et aLl 
I2012at tHosokawa et al. 201 1; Clark et al. 2(3lB. It is even 
likely that the Pop III IM F can be considerably va riable 
depending on environment dO'Shea & Normanll2007h . We 
choose not to explore various analytic forms for the IMF, as 
constraining it is beyond the scope of this paper Instead, we 
consider that the main way by which the Pop III IMF can in- 
fluence galaxy formation, in contrast to the known Pop II IMF, 
is by enhancing the output of ionizing radiation and the num- 
ber and intensity of supernovae. In particular, PISNe, which 
are hypothetically plausible from stars in the mass range 
140-260 M(7) (or for lowe r masses if rotation is considered, 
see e.g. IStacv et al.l 12012b). would potentiall y be dramatic 
singular ev ents in the evo lution of any galaxy (I Bromm et alj 
2003; Whal en et al.ll2008l) . To account for the occurrence of 
PISNe we use two different particle masses for Pop III stars. 
Every newly formed Pop III stellar particle is randomly as- 
signed to be either a I7OM0 star, which is to explode in a 
PISN, or IOOM0 star, which only generat es a mild exp losion 
before collapsing into a black hole (Heg er & Woosley||2d02i 
I2OIOI) . The proportion of these two types of particle mass and 
fate is governed by a single parameter, PpisN, which is the 
fraction of PISNe progenitors (170Mo stars) that form when 
the Pop III star formation criteria are met. In our fiducial runs, 
we set PpjsN = 0.5. This value was chosen to test the maximum 
possible impact of PISNe on galaxy evolution, and probably 
represents the most top-heavy the primordial IMF can pos- 
sibly be. Since the atmospheres of Pop III stars are free of 
metals, they are unable to drive stellar winds and therefore do 
not enrich the ISM in any way other than through supernovae. 
Pop III stars which have masses too low to produce SNe are 
ignored in our model. 

For our fiducial value of nn.mm (as well as all other param- 
eters considered in Section |2.3.1t , we found that the gas mass 
in a maximally refined cell at z ss 20 is sometimes not suf- 
ficient to form a I7OM0 star. Therefore, we prevent further 
refinement in metal-free cells that have nn > 0.5nH,mm and 
whose splitting would leave insufficient mass to form the star 
Through tests, we have checked that this refinement restric- 
tion never artificially slows down Pop III star formation. It 
becomes especial ly rele vant in the super-Lagrangian runs dis- 
cussed in Section l2.3.3l and in the H and Q Mpc boxes which 
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inherently have very high resolution. 

The PISN from a HOM© star releases 27 x 10^' erg of 
thermal energy, as well as 80 M© of metals into th e ISM 
(iHeger & Woosley 2002). As suggested by Wise & Abell 
(1^08), we use a delay of 2.3 Myr from the formation of a 
170Mq particle to its PISN event, representing th e main se- 
quence lifetime of this type of star (i S chaereiil2002l) . After the 
supernova goes off, the cell which hosts it often winds up with 
super-solar metallicity. This presents problems for the cool- 
ing functions employed in our code, so we turn off all metal 
cooling for gas at temperatures higher than 10* K. Accord- 
ing to the models of Heger & Wooslev (2002), a 170Mq star 
is completely disrupted by its PISN event, and all gas mass 
from the stellar interior would be ejected into the ISM leaving 
no remnant. The ejecta then consists of 80 M© of metals from 
the core, as well as the primordial envelope which is 22 M© 
of He and 68 Mq ofH. 

A IOOMq star does not explode as a PISN, but its actual 
fate is still uncertain and depends on the details of the stel - 
lar rotation and magnetic structure (iHeger & Woosleyll2010^ . 
Before undergoing core collapse, such a star would experi- 
ence thermonuclear pair-instability pulses that would eject the 
outer layers of H and He, with possible traces of the elements 
C, N, and O. The energy released in such pulses is of the or- 
der or smaller than the energy of a normal SN type II. If the 
remaining core has enough rotation to trigger a ga mma-ray 
burst in the collapsar model (IWooslev & HegeJl2012l) . it may 
then lead to a powerful explosion with over 10^2 erg of en- 
ergy and the ejection of significant mass of metals. Without 
rotation, the core may drive a weak collapsar explosion or 
no explosion at all, when all the remaining mass recollapses. 
Given these uncertainties, and to contrast with the case of 
a full PISN explosion for the 170Mq stars, for the IOOMq 
stars we assume that no substantial metals are deposited into 
the ISM and that the released energy corresponds to a stan- 
dard SN type II. About 50 M© of gas is released into the ISM, 
while also leaving behind a remnant black hole of -50 Mq. 

To prevent artificial radiative losses, PISN energy and mass 
ejecta are distributed within a sphere of constant density, with 
a radius 1.5 cell lengths, centered at the middle of the PISN 
host cell. Each of the 27 cells within such a sphere, consist- 
ing of the star's host cell and its immediate neighbors, receive 
a dose of energy and metal-rich gas proportional to the vol- 
ume of the cell contained within the sphere. This prescription 
is physically motivated as we found that our typical timestep 
(about 450 yr) is too coarse compared to the typical timescale 
of the early free expansion phase of the SN remnant. For ex- 
ample, it takes the ejecta -500 yr to traverse half of the typical 
Pop III star host cell length (4.5 pc) at z-20 if it travels at the 
free-expansion velocity. This velocity is computed here by 
assuming that all of the PISN energy of 27 x 10^' erg goes 
into kinetic energy of the ejecta. 

In addition to the supernova feedback, all Pop III stars have 
enhanced radiative feedback relative to Pop II counterparts, 
due to the lack of metals in their atmospheres ( Schaerer 2002). 
We enhance the radiative output of Pop III stars by a factor of 
10 relative to Pop II, following Wise & Cen (2009). However, 
after a Pop III stellar particle undergoes supernova, radiative 
feedback from the star is completely shut off. On the other 
hand. Pop II stellar particles shine according to a light cu rve fit 
from Starburst 99 model results (Leithe rer et al.|[l999l) . This 
light curve consists of a flat component for the first 3x10^ 
yr, followed by a steep power-law falloff. Radiative feed- 



back from Pop II stellar particles becomes insignificant after 
3x10^ yr. Since Pop II stars shine longer than both types of 
Pop III stars, the factor of 10 radiative enhancement does not 
translate into a proportionally higher number of ionizing pho- 
tons per lifetime. Pop II stars emit 6,600 ionizing photons per 
stellar baryon per lifetime. In our fiducial runs. Pop III stars 
emit 38,800 and 34,500 photons per baryon per lifetime for 
the lOOM© and 170Mq stars, respectively. 

Pop III stars form in gas that has metallicity logjyZ/Z© < 
-3.5. This threshold is chosen to match the critical metallicity 
discovered by Bromm et al. (2001), and has held up in later 
studies (ISmith et al.l 12009). Though the exact value of this 
critical metallicit y is still uncertain and can be affected by the 
presence of dust dOmukai et al.ll20()5h . we find that it is not 
very important as the majority of Pop III stars form in truly 
primordial, or nearly primordial gas. Compiling all of our 
simulations, we found that only -10% of Pop III stars form 
with logjoZ/Z© > -5. 

We summarize the formation criteria for Pop III stars with 
the following set of equations, 

«// > "H./mn = 10"* Cm"^ 
/h2 > .fH2.min = 10 

logioZ/Zo<-3.5. (3) 
Values given for each variable represent the fiducial choices. 

2.3. Convergence Study & Setting Fiducial Parameters 

In this section, we describe the test runs that justify the nu- 
merical setup and the choice of parameters for our main runs. 
Since the Pop III star formation recipe described above is one 
of the critical components of our study, we focus on testing 
the key elements of this model. In Table[T]we list the details 
of the simulations performed in our suite. Box OverDense has 
many more potential sites for Pop III star formation than the 
other 1 /;"' Mpc boxes, and therefore serves as the best test- 
ing ground. It is important to keep in mind that while every 
parameter we test has an effect on Pop III star formation, the 
most drastic differences between the simulations are caused 
by the choice of initial conditions. The role of cosmic vari- 
ance will be explored more comprehensively in Paper II. 

2.3.1. Density threshold for Pop III star formation 

First, we choose an appropriate value for the density thresh- 
old for creating Pop III stars. In Figure [T] we examine the 
high-density end of the volumetric probability distribution 
function (PDF) of the hydrogen number density at a = 0.085 
(z = 10.8). In order to test the properties of the primor- 
dial gas from which the first stars form, we ran a special 
set of simulations with no star formation or chemical en- 
richment (runs OverDense_noSF_8, UnderDenseH-_noSF_8, 
UnderDense-_noSF_8, as well as additional versions of each 
with one more and one fewer maximum level of spatial re- 
finement). Though the total mass of gas in each box is the 
same, only gas in the most massive halos has collapsed to 
this density regime, meaning that the PDFs are very sensitive 
to the number and nature of such halos. The PDFs of Box 
OverDense and UnderDense-n are offset by a factor of -5 at 
all densities, while the UnderDense- box is offset from Box 
UnderDense-H by another factor of -5. This difference is also 
seen in the maximum density achieved in each box. For our 
fiducial resolution of 8h-8 levels, all three boxes are able to 
reach a density of at least 10000 cm""* by a = 0.085, giving us 
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Table 1 

Simulation runs 



Run Base grid imax dx (pc) nioA/CM©) nH,mm{cm^) Description 



Convergence study 



UnderDense-_noSF_7 


2563 


7 


44 


5500 


No star formation 


UnderDense-_noSF_8 


2563 


8 


22 


5500 


No star formation 


UnderDense-_noSF_9 


2563 


9 


11 


5500 


No star formation 


UnderDense+_noSF_7 


2563 


7 


44 


5500 


No star formation 


UnderDense+_noSF_8 


2563 


8 


22 


5500 


No star formation 


UnderDense+_noSF_9 


2563 


9 


11 


5500 


No star formation 


OverDense_noSF_7 


2563 


7 


44 


5500 


No star formation 


OverDense_noSF_8 


2563 


8 


22 


5500 


No star formation 


OverDense_noSF_9 


2563 


9 


11 


5500 


No star formation 


OverDense_noSF_HMpc 


2563 


8 


11 


690 


No star formation, 0.5 ft 



Fiducial runs 



UnderDense-_nHle4_fid 2563 8 22 5500 10000 Underdense box with no massive galaxies, fiducial parameters 

UnderDense+_nHle4_fid 2563 8 22 5500 10000 Underdense box with one massive galaxy, fiducial parameters 

OverDense_nHle4_fid 2563 8 22 5500 10000 Overdense box, fiducial pai'ameters 

Density threshold study 



OverDense_nHle3 2563 8 22 5500 1000 Lowest density threshold for Pop 111 star formation 

OverDense_nH5e3 2563 g 22 5500 5000 Low density threshold for Pop 111 star formation 

OverDense_nH2e4 2563 8 22 5500 20000 High density threshold for Pop III star formation 



Mass resolution & refinement criteria 



OverDense_SL7 


2563 


8 


22 


5500 


10000 


Super-Lagrangian refinement 0.7* 


OverDense_SL5 


2563 


8 


22 


5500 


10000 


Aggressive super-Lagrangian refinement 0.5* 


OverDense_HiRes 


5123 


7 


22 


690 


10000 


Higher mass resolution 


OverDense_HMpc 


2563 


7 


22 


690 


10000 


0.5 ft"' Mpc box 


OverDense_HMpc_HiRes 


2563 


8 


11 


690 


10000 


0.5 ft"' Mpc box, higher spatial resolution 


OverDense_HMpc_SL5 


2563 


7 


22 


690 


10000 


0.5 ft"' Mpc box, super-Lagrangian refinement 0.5* 


OverDense_QMpc 


2563 


8 


5.5 


86 


10000 


0.25 ft-' Mpc box 



Alternative physics 



OverDense_ExtremeSN 2563 8 22 5500 10000 Extreme PISNe (Section lTT6l 

OverDense_ExtremeRad 2563 8 22 5500 10000 Extreme Pop 111 radiation field (Section [TX6l 

OverDense_LowMass 2563 8 22 5500 10000 Pop III IMF and feedback mirror Pop II (Section|233) 



Column 1.) Name of the run; 2.) Base grid, number of DM particles, number of root cells; 3.) Maximum number of additional levels of refinement; 4.) Minimum 
cell size at the highest level of refinement in comoving pc; 5.) DM particle mass in Mq; 6.) Minimum H number density for Pop 111 stai' formation in cm-3; 7.) 
Further description of the run. 



enough time to study star formation in every box before our 
stopping point of a = 0. 1 . 

Based on these results, we chose «H,mm = lOOOOcm"^ as our 
fiducial value for the density threshold. In addition to the con- 
straints obtained from the PDF, other considerations went into 
this selection. A lower value would suffice to meet our proto- 
cloud collapse criteria, but would result in all Pop III stellar 
particles forming before cells are maximally refined. Such 
an outcome is poor practice in hydrodynamic simulations, as 
subgrid physics is being invoked on scales where the resolu- 
tion is still good enough to self-consistently capture relevant 
physical processes. On the other hand, using a higher thresh- 
old would allow the maximally refined cells to reach densities 



beyond the resolving power of the simulation. When such 
conditions are reached, either further refinement or subgrid 
physics should already be in use. In addition, using a higher 
density threshold in our test runs often led to Pop III stars 
forming in bursts (in the same timestep, in neighboring cells). 
We do not speculate here whether such bursts are physically 
plausible or not, but the scales necessary to model this process 
properly are certainly unresolved in our simulations, and we 
suspect that higher temporal or spatial resolution would re- 
veal that feedback from the first star in a cell would suppress 
further clustered star formation. 

To determine the ultimate effect of the density threshold on 
Pop III star formation, additional runs were performed with 
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Figure 1. The distribution function of hydrogen number density for runs 
without star formation at a = 0.085 (z = 10.8). Blue lines are for Box 
UnderDense-, red for Box UnderDense+, black for Box OverDense. Dotted 
lines are for 8+7 levels of refinement, solid lines for 8+8, short-dashed lines 
for 8+9. The long-dashed green line represents our chosen density thresh- 
old ^Hjiiin — 10000 cm ^ . All runs with at least 8 levels of refinement have 
sufficiently dense gas to form stars by this epoch. 



nHjnin = 1000, 5000, and 20000 cm"^ using the Box Over- 
Dense initial conditions. Varying this threshold by a factor 
of 20 changes the scale factor at which the first star forms 
only from a = 0.0463 to a = 0.0483, or from redshift z = 20.6 
to z = 19.7. In the n/f_„„„ = 1000 cm""* run, the density thresh- 
old is reached at a lower level of refinement from the other 
cases considered, allowing the first star to form sooner The 
variation in the other three runs is only from a = 0.0480 to 
a = 0.0483. Such marginal differences demonstrate that for 
a given set of initial conditions, our actual density threshold 
criterion for the formation of Pop III stars has little effect on 
when and where they form. In each of these runs, only a single 
Pop III star formed in each box before a = 0.05, and the total 
number of Pop III stars varied between 4 and 5 at a = 0.055. 
Based on these tests, we have determined that the total number 
of Pop III stars formed had little correlation with the density 
threshold within the considered range. 

After the first star forms in a given halo, the gas den- 
sity can be significantly reduced near the center, quenching 
further star formation. Figure |2] demonstrates this effect in 
run OverDense_nHle4_fid. The PDF of this galaxy is de- 
pleted at high density 10 Myr after a PISN explosion. The 
corresponding galaxy from the run without star formation, 
OverDense_noSF_8, is also shown for reference. While the 
galaxy in OverDense_noSF_8 contains some dense gas above 
riH = 10 cm"^, it is depleted in our fiducial run. Since this den- 
sity is nowhere near any nn.mm that we have considered in our 
tests, we can conclude that Pop III stars will not form in quick 
succession in this halo. 



2.3.2. Molecular Fraction 
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Figure 2. PDF for most massive galaxy in Box OverDense at a = 0.05 (z = 
19) with star formation (solid, ran OverDense_nHle4_fid) and without star 
formation (dotted, run OverDense_noSF_8). The Pop III star that recently 
formed at the center of this galaxy has temporarily depleted it of the dense 
gas necessary for more star formation. 



Another component of the Pop III star formation criterion 
is the requirement of a minimal fraction of H2 in the host cell. 
To determine what value of the H2 threshold makes sense in 
the context of these simulations, we examine the molecular 
fraction of hydrogen as a function of density in the runs with- 
out star formation at the epoch (z w 20) when gas is begin- 
ning to reach densities close to rtH.min- Figure [3] shows that 
the molecular fraction in primordial gas generally increases 
with density, but saturates above «// ~ lOcm"^. The satura- 
tion value of ///, grows slowly with time in the absence of 
star formation, and does not appear to depend significantly on 
spatial or mass resolution. Our fiducial choice of 10"^ for the 
minimal H2 fraction does not exclude dense gas from forming 
stars in any runs, as long as little Lyman- Werner radiation is 
present. 

2.3.3. Super- Lagrangian Refinement 

Since we have demonstrated resolution dependence for the 
maximum density of gas within a given galaxy, it is expected 
that refinement criteria could play a role in controlling when 
gas in the simulation first reaches the nH.m/n threshold. To test 
this, in some of our runs we employ super-Lagrangian (SL) re- 
finement criteria. The refinement criteria in a cell can be writ- 
ten as 2 X X x x 0.8 for the dark matter mass, and 

0.3 X niDM X xX^ X 0.8 for the gas mass, where £ is the 
level of the cell which is to be refined. In this formalism X = l 
implies Lagrangian refinement as described at the beginning 
of Section |2] We have tried runs with very aggressive SL re- 
finement (X = 0.5, run OverDense_SL5) and less aggressive 
refinement (X = 0.7, run OverDense_SL7). Running these 
simulations in Box OverDense with nH.min = lOOOOcm"-^, we 
found that the epoch at which Pop III stars first appear is 
pushed back from a = 0.0478 to a = 0.0456 with X = 0.7, and 
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Figure 3. Molecular fraction of hydrogen vs number density for three 
Box OverDense runs using different resolution without star formation at 
a = 0.05 (z = 19, runs OverDense_noSF_9, OverDense_noSF_8, and Over- 
Dense_noSF_7 are black, red, and blue, respectively). The median H2 frac- 
tion in each density bin is indicated by a triangle, while the error bars show 
25th and 75th percentile levels. At this epoch, when the first stars would nor- 
mally be forming, our fiducial resolution of 8 levels has the same molecular 
fraction as if we were using one more or one fewer level of refinement, and 
the points actually lie directly on top of each other for rifj < 10 cm"''. Run 
OverDense_noSF_HMpc at a = 0.05 (green) has lower values and a relatively 
wider spread of H2 fraction for iifj < crrT^ but converges with the other 
runs at higher densities, demonstrating a lack of dependence on mass reso- 
lution. Also shown is run OverDense_nHle4_fid where star formation has 
already taken place by a = 0.05 (purple). Since the gas has been enriched 
by a PISN, molecular gas can form at much lower densities. Light red and 
light blue points trace out the H2 fraction in every single cell for run Over- 
Dense_noSF_8 and OverDense_nHle4_fid, respectively. 



to a = 0.0435 with X = 0.5. This demonstrates that the use 
of SL refinement is an important numerical tool for explor- 
ing the earliest epoch of star formation in a given simulation 
box. However, using SL refinement produces an enormous 
number of high-level cells: at a = 0.05, run OverDense_SL5 
has a factor of 4000 more maximally refined cells than run 
OverDense_nHle4_fid. This drastic difference makes the SL 
simulations prohibitively expensive soon after the first stars 
form. 

Therefore, we use these SL runs to study Pop III star for- 
mation at the earliest possible epochs, when the mass of the 
halos that hosted them was low enough for PISNe to have 
their maximal effect. 

2.3.4. Increased Mass Resolution 

We test the effects of mass resolution by setting up one run 
with 512^ initial grid, giving a DM particle mass of 690 M©. 
We use 7 additional levels of refinement, therefore granting us 
the same maximum spatial resolution as in the fiducial 256^ 
run. Having consistency in spatial resolution allows us to test 
the effects of mass resolution alone. All other numerical pa- 
rameters are kept consistent with run OverDense_nHle4_fid. 

The increased mass resolution has several immediate im- 
plications. Since we now resolve halos of mass IQ^Mq with 



over 1000 particles, we can better probe the regime where 
the very first Pop III stars are expected to collapse in proto- 
galactic minihalos. Indeed, the epoch of formation of the first 
star in the box becomes a = 0.0427 in a halo of L5 x IO^'Mq 
(compared to a = 0.0481 and M/, = 7.5 x 10'' in run Over- 
Dense_nHle4_fid). The higher mass resolution effectively 
means that there is more power on the small scales respon- 
sible for the growth of halos in this mass regime. Due to the 
high computational cost of this run, we have only advanced it 
too = 0.055. 

The effect of increasing mass resolution is further explored 
through runs using the H and Q boxes of 0.5 Mpc and 0.25 
/i"' Mpc in size. Applying the 256"' base grid to these boxes 
gives us a DM particle mass of 690 M0 and 86 Mq, respec- 
tively. We find that the H box (run OverDense_HMpc_HiRes) 
produces a Pop III star by a = 0.0456 in a halo of mass 
1.5 X 1O*'M0, while the Q Mpc box (run OverDense_QMpc) 
does not produce one until a = 0.0473. Using SL refinement 
in the HMpc box (run OverDense_HMpc_SL) allows us to 
see a Pop III star forming in a 8 x 1O^M0 halo. 

The earlier formation epochs and lower mass of halos host- 
ing the first stars in the H and Q Mpc boxes, compared to the 
fiducial 1 /;"' Mpc runs, show that it is crucial to have high 
enough mass resolution to capture Pop III star formation in 
halos close to 1O''M0. It has been previously shown that ha- 
los less massive than this threshold will not achieve significant 
enough H2 abundances to t rigger Pop III star formation at ear- 
lier epochs ( lYoshida et al.ir2003 ). The further significance of 
this mass range will be explained in our Results section. Fig- 
urelHshows explicitly how varying refinement criteria, spatial 
resolution, and mass resolution affected the lowest possible 
mass for a star-forming galaxy. 

2.3.5. Low Mass Pop 111 IMF 

We present one simulation, run OverDense_LowMass, 
which does not rely on a top-heavy IMF for Pop III stars. The 
conditions for Pop III star formation in this run are similar to 
our fiducial top-heavy recipe in that we use the same thresh- 
old rtH./nin to determine which cells are allowed to form stars. 
However, the stellar particle masses are drawn from the same 
IMF as for Pop II stars. This run explores the possibility that 
Pop III stars were ordinary low-mass objects. Whenever the 
density in a given cell exceeds the threshold, the star forma- 
tion rate is determined according to the following relation: 



P* ~ Pgas, 
TSF 



(4) 



where pgas is the mass density of all gas in the cell. This 
relation is similar to equation IH but does not explicitly use 
molecular hydrogen. This modification is necessary because 
primordial gas can reach densities above our star formation 
threshold, but cannot become fully molecular without the 
presence of dust. 

2.3.6. Extreme Pop 111 Feedback 

To isolate the relative impacts of the feedback effects, we 
ran toy simulations using exaggerated values for the PISN 
energy and ionizing photon yield. In one run, called Over- 
Dense_ExtremeSN, PISNe released 270 x 10^' erg of thermal 
energy, a factor of 1 larger than in all other runs . The extreme 
ionizing simulation OverDense_ExtremeRad had instead an 
additional factor of 10 boost in the ionizing photon flux of 
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Figure 4. The least massive galaxy to host a star in various runs vs. the 
minimum comoving cell size employed in the run. The colors indicate the 
simulation box size: 1 /i-' Mpc (black), 0.5 ft"' Mpc (blue), or 0.25 /i"' Mpc 
(red). The shape of points indicates the refinement criterion employed in the 
box, with open squares for Lagrangian refinement, filled triangles for 0.7^ SL 
refinement, and open five-pointed stars for 0.5^ SL refinement. 



Pop III Stars, giving the 100 M© and ITOMq stars 388,000 
and 345,000 photons per lifetime, respectively. While these 
models are too strong to be consistent with any published re- 
sults, using them allows us to explore the most extreme effects 
of Pop 111 feedback. 

3. RESULTS 

Pop 111 stars in our simulations begin to form in halos of 
mass Mh > lO^'M© starting at a « 0.045 (z ft! 21.2) in ac- 
corda nce with expectations from prior work dYoshida et alJ 
l2003h . Figure |5] shows the mass of host halo in which each 
Pop 111 star formed in our fiducial runs for the three 1 /i"' Mpc 
boxes and the H Mpc box. In the 1 /i"' Mpc runs, the halo 
mass for first star formation is close to lO^M©. This mass 
is an order of magnitude larger than that of the halos h ost- 
ing the firs t stars in the simulations of I Wise et alj (12012 1) and 
iGreif et d] (1201 Ih, and those p referred by theoretical consid- 
erations dTegmark et alj|1997[) . Consequently, those authors 
also find an earlier epoch for the formation of the first stars. 
Given that the extra mass resolution granted by the H Mpc box 
allows us to see star formation in 10^ Mq halos, we acknowl- 
edge that our fiducial 1 Mpc runs are not properly re- 
solving the very first star-forming minihalos. Rather, they are 
more generally simulating Pop 111 star formation in an early 
population of galaxies. The fraction of star-forming halos in 
run OverDense_nHle4_fid at a = 0.07 (z= 13.3) was only 1% 
in the mass range 10* Mq <M/, < IO^Mq, but it reached 65% 
for Mh > lO^M©. In run OverDense_HMpc_HiRes, these 
numbers increase significantly to 15% for 10*' M© < M/, < 
lO^M© and 100% forM;, > lO^M©. 

It is also worth noting that the ratio of star-forming galax- 
ies in the range lO^M© < M;, < lO^M© falls off gradually 
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Figure 5. Each Pop III star's host halo mass at the time of formation vs the 
scale factor at which the star formed for run UnderDense-_nHle4_fid (blue), 
run UnderDense+_nHle4_fid (red), run OverDense_nHle4_fid (black), and 
run OverDense_HMpc_HiRes, which did not go past a = 0.075 (green). Ad- 
ditional points for the halos hosting the first stars from each of the runs used 
for Figure |4] are also included, with the same color and shape scheme. In 
our 1 /j"' Mpc runs, Pop III star formation happens almost exclusively be- 
tween 10' Mq and 10* Mq. The additional mass resolution in run Over- 
Dense_HMpc_HiRes makes it possible to see that the first Pop III stars really 
form in halos between 10* Mq and IO'Mq. The average mass of Pop III 
star-forming halos increases slightly with time. When multiple Pop III stars 
form within a galaxy in a very short time interval, points on the plot are 
grouped into a clustered shape. 



with time in run OverDense_nHle4_fid. It changes from 
65% at a = 0.065 to 20% at a = 0.1, suggesting that halo 
mass alone is not a good proxy for determining whether a 
galaxy can achieve the high density required for our Pop 111 
star formation criteria. One potential cause of the change is 
the decreased physical spatial resolution at later epochs, but 
according to our study of the gas in the first star-forming 
galaxy shown in Figure [51 the factor-of-two difference in 
spatial resolution granted by using one fewer level of re- 
finement does not preclude gas from reaching the fiducial 
"H./nm = 10000 cm"^ threshold. The difference is more likely 
to be rooted in the variation of physical density in halos of a 
given mass with time. For halos between 10^ and 10** M©, 
the average matter density within the virial radius changes 
from 2.5 x lO'^Mepc'^ at a = 0.065 to 6.8 x 10-^M©pc-'' 
at a = 0.1, due to the expansion of the universe. Even the 
density within the central 100 pc of these halos changes from 
0.73 M© pc""* to 0.34M© pc"-' between the same two epochs. 

Very few Pop 111 stars formed in halos with Mi, > 10**M©, 
because such halos have akeady been enriched to metallic- 
ities above logjgZ/Z© = -3.5, allowing for normal star for- 
mation to kick in. One major exception occurs in run Over- 
Dense_nHle4_fid in a galaxy that has already formed a sig- 
nificant number of Pop 11 stars that have in turn enriched the 
ISM terminating further Pop 111 star formation. However at 
a = 0.095, this galaxy undergoes a major merger with another 
massive halo, causing low-metalHcity gas in the outer part of 
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Figure 6. Though metals from PISNe are ejected past the virial radius, they 
do not stay there for long. Plotted here are the radii enclosing 80% of the 
metals produced in the galaxy, as well as the radii where the enclosed mass 
of baryons divided by the virial mass equals 80% of the universal baryon 
fraction. Each line represents a galaxy as it evolves in time, beginning at the 
epoch when the first star forms. Within 50 Myr of the PISN event (denoted 
by five-pointed stars), most of the gas and metals have begun to recollapse, 
or are at least enclosed within the virial radius once again. The maximum ex- 
tent of metal propagation is strongly regulated by galaxy mass. The distance 
between two points on each line corresponds to 10 Myr. 
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Figure 7. The evolution of the baryon fraction vs. the evolving halo mass. 
Each line traces a galaxy from the time of formation of the first star. The 
baryon fraction is computed within the virial radius, and is normalized by 
the universal value /j, ,„„. While supemovae initially cause a depletion of 
baryons in galaxies of M/, < 3 X 10' Mq, this depletion is only temporary. 
In galaxies of Mj, > 3 X lO' Mq, there is no strong evidence that PISNe are 
able to deplete baryon fractions. Notation is the same as in Figure|6] with the 
addition of triangles to represent SN II produced by IOOMq Pop III stars. 



the galaxy to collapse, thereby triggering a burst of Pop III 
star formation. All of these stars form with metallicities right 
around the critical logjgZ/Z© = -3.5 value, suggesting that 
their existence is sensitive to the value of this threshold and 
therefore should not be treated as a general result. 

3.1. Effect of Pop III stars on their host galaxies 

The strong ionizing flux of Pop III stars and enormous en- 
ergy injections from PISNe have been shown in previous work 
to significantly alter the ISM of their host galaxies. Here we 
explore such effects during time when Pop III stars are the 
dominant drivers of feedback. 

In Figure |6] we show that there is a significant variation in 
the potential effect of Pop III stars on their host galaxies de- 
pending on the galaxy mass. In halos with M/, < 3 x IO^Mq, 
Pop III stars can temporarily evacuate the gas from the galaxy, 
and the metals from PISNe are ejected past the virial radius 
into the intergalactic medium (IGM). On the other hand in ha- 
los with Mh > 3 X 1O''M0, the metals are confined within the 
virial radius, and there is little movement of baryons beyond 
the virial radius. 

This divide can also be seen in Figure |7] which shows the 
baryon fraction computed within the virial radius for thirteen 
halos taken from a variety of runs. Again, a noticeable thresh- 
old at Mh = 3 X 10'' M0 distinguishes galaxies that have their 
gas evacuated by PISNe from those that do not. Galaxies less 
massive than this threshold typically have their gas content 
plummet by at least a factor of two within 10-30 Myr after the 
PISN, with the least massive ones falling below 10% of the 
universal baryon fraction. In contrast, more massive galax- 



ies lose a much smaller percentage of their gas and end up 
with baryon fractions in excess of their pre-explosion values 
within -100 Myr This dividing line between "low-mass" and 
"high-mass" galaxies is therefore a logical choice for param- 
eter that distinguishes different regimes of Pop III feedback. 
We examine these two regimes separately below. 

3.1.1. Dependence on halo mass 

To further understand how halo mass can determine the ef- 
fectiveness of Pop III stellar feedback, we examine the struc- 
tural evolution of several galaxies in different mass regimes. 

First, we examine a relatively low-mass galaxy from run 
OverDense_HMpc_HiRes, which is shown by the black line 
in Figure|6l and the black fine that extends to M = 9 x 10^ Mq 
in Figure]?] The first star (of I7OM0) forms when the mass of 
the halo is 1.7 x 1O''M0. Within 20 Myr of its formation, the 
PISN has blasted metals out beyond 1 kpc from the galactic 
center (or 4/?,,,> at this epoch), and the baryon fraction has 
dipped to fb/ fb.uni ~ 0.15. However, soon after this point the 
baryon fraction begins to grow again, and the virial radius 
increases enough to enclose a larger fraction of the expelled 
gas and metals. 

At f = 1 39 Myr after the first PISN (the halo mass is now 9 x 
IO^'Mq), the galaxy has regained -37% of the PISNe metals. 
The baryon fraction is over half of the universal value. It will 
still take more time for this galaxy to completely recover from 
the explosion, but there is considerable evidence from Figure 
]6]that metals and baryons in general are flowing in rather than 
out of the galaxy. Another sign of recovery is that Pop II 
star formation has commenced within the galaxy, as it now 
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Figure 8. The ejecta of PISNe is traced via examining tlie enclosed mass 
of metals as a function of galactocentiic distance. Lines of different color 
correspond to 1, 27, 59, 91, and 139 Myr after the first PISN. A second PISN 
happens 15 Myr after the first in a nearby galaxy. Approximately 60 Myr after 
the first PISN, more metals are flowing into the galaxy than outwards, as the 
metal-rich ejecta have mixed with primordial gas accreting onto the galaxy. 
Arrows show the direction of metal movement at each epoch. The length of 
each an'ow corresponds to the distance traversed by the metals in a 20 Myr 
interval. The y-axis positions of the arrows show the mass of metals at each 
epoch used to compute the rate of propagation. This galaxy is depicted by the 
black line in Figures |6]and|3 



contains 4 Pop II stellar particles (which still contribute little 
to the metal budget). 

Figure [8] follows the radial distribution of metals in this 
galaxy from the time of the first PISN until 139 Myr after 
it has exploded. The Pop II stars have contributed less than 
1 M0 to the metal budget, so essentially all of the metals 
shown here are products of the first PISN, and of a second 
PISN which happens 15 Myr after the first in a neighboring 
halo at a distance of approximately 2 kpc. The mass of this 
galaxy has increased by a factor of ~5 between a = 0.0508 
when the star first formed and a = 0.0726 at the final snapshot 
considered. Such rapid growth bring the metals back into the 
galaxy. 

While the PISNe do clearly cause baryon depletion and sup- 
press star formation in galaxies, the rate of growth of these 
galaxies is large enough that a mixture of ejecta and fresh pri- 
mordial gas from the IGM fall in to restore the baryon frac- 
tion to at least 50% of the universal fraction within -150 Myr 
This replenishment results from a combination of actual re- 
accretion of ejected material as seen in Figure[8] and the rapid 
growth of galaxy virial radius. Ultimately, Pop III stars have 
little direct impact on the IGM. 

Next, we study a galaxy from run OverDense_nHle4_fid 
that was the first to form a Pop III star, which also happens 
to be a PISN progenitor. This galaxy is represented by the 
brown lines in Figures |6] and [7] The first star forms when 
the mass of the halo is 7.5 x lO^'M©. About 35 Myr later, 
80% of the metals generated in the PISN have propagated as 
far as 525 parsecs from the core. The injection of metals by 



the PISN is enough to bring the gas metallicity in some cells 
hundreds of parsecs away from the galactic center to be as 
high as logjgZ/Z© = -1. After another 50 Myr, the effects 
of the outflow have subdued. Not only are 80% of generated 
metals now entirely confined to the innermost 120 parsecs, 
but the metals have diffused, and the maximum metallicity 
has decreased by 1 dex. This suggests that the inflow of new 
primordial gas is playing a greater role in the evolution of the 
galaxy than the outflow generated by the PISN. The majority 
of metals produced by PISN do not escape into intergalactic 
space. 

We emphasize that some galaxies in this mass range should 
have hosted Pop III stars at earlier times than those resolved 
by our simulations. The apparent ineffectiveness of Pop III 
feedback demonstrated here shows that simulations which do 
not resolve halos with M/, < 3 x lO^M© are missing a por- 
tion of galactic evolution. This omission could mean that Pop 
III stars should self-terminate at earlier times, hence decreas- 
ing their contribution to the cosmic ionizing background. On 
the other hand, the expulsion of baryons from low-mass halos 
leads to suppression of Pop II star formation, which implies 
we may also overestimate the Pop II rates. The balance of 
these effects will be explored further in Paper II. 

Figure |9] shows how far metals propagate in galaxies rel- 
ative to the stellar cores. The "gas metal half-mass radius" 
is calculated as the more familiar stellar half-mass radius, but 
tracing the total mass of metals in the gas phase instead of stel- 
lar mass. In general, the metals are always able to propagate 
out beyond the stellar cores, but the extent depends strongly 
on galaxy mass. For a homogeneous comparison, and po- 
tential future probes by observation, we plot all star-forming 
galaxies from a single epoch, a = 0.075 (z = 12.3). In galax- 
ies with Ml, > 1O^M0, the metals remain within a factor of 
2 of the stellar radius. In less massive galaxies, metals are 
able to propagate further, sometimes by as much as a factor 
of 10, owing to the lower potential wells of these galaxies. 
Nonetheless, considering that the stellar half-mass radii of all 
our galaxies range 5-30 pc, the location of the bulk of metals 
is still limited to only the innermost regions of galaxies. 

3.2. Effects of the uncertainty in Pop III feedback and IMF 

With our additional runs, we can check if the relatively 
inefficient feedback is due to the specific fiducial parame- 
ters that we adopted. However, even with th e extre me Pop 
III feedback prescriptions described in Section l2.3.6i we find 
that the baryon fraction within the virial radius is never sig- 
nificantly depleted. At a = 0.05 (z = 19), 10 Myr after the 
PISN explosion in the first star-forming galaxy, we find ft, = 
9.6% and 10.7% in runs OverDense_ExtremeSN and Over- 
Dense_ExtremeRad, while it is 11.5% in the fiducial run. This 
is a relatively small resulting difference for a 10-fold increase 
in the thermal energy and ionizing radiation output. The mass 
of the galaxy at this epoch is 9 x IO^Mq, which we have 
shown to be large enough to withstand standard Pop III feed- 
back. 

The effect of extreme feedback on the propagation of metals 
is stronger Metals tend be blown out of galaxies anisotrop- 
ically, often extending outwards in directions orthogonal to 
filaments, into lower density regions. The galactocentric ra- 
dius that encloses 80% of the metals formed in the PISN 
stretches out to 2.0/?i,,> 20 Myr after the explosion in run Over- 
Dense_ExtremeSN, while this radius reaches just 1.3/?,,,> in 
the fiducial run. These metals do not fully escape the gravita- 
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Figure 9. The ratio of gas metal half-mass radius to stellar half-mass radius 
vs halo mass for galaxies which have had at least one PISN, at a = 0.075 (z = 
12.3) in runs OverDense_nHIe4_fid (black) and OverDense_HMpc_HiRes 
(green). Metals propagate further relative to the stellar cores in galaxies of 
lower mass. 



tional pull of the galaxy, however, and 90 Myr after the explo- 
sion, the galaxies from both extreme feedback runs contain 
80% of the metals from the first PISN within the virial radius 
(in the fiducial run, they are contained within just 0.15/?,,,/)- 
The PISN metal ejecta in OverDense_ExtremeRad followed a 
roughly similar trajectory to those in OverDense_nHle4_fid. 

Though we increased the feedback effects by a factor of 
10, only modest and temporary differences were observed be- 
tween the runs. Such inefficiency of feedback demonstrates 
the weak coupling of thermal energy from PlSNe to the ISM 
at the densities and temperatures resolved by our simulations, 
as almost any amount of energy can be quickly radiated away. 
This can be seen when considering typical cooling times in 
the ISM, Tcooi = khT/A n fti 3000(^/10^ K)(l cm-Vn)yr, for 
A = 10"^^ ergs"' cm"^ (iHopkins et alJ |20TT1) . This cooling 
time is comparable to just a few typical timesteps in our sim- 
ulations for cells at these densities (-500 yr). The impacts 
of extreme feedback were more pronounced in the IGM, par- 
ticularly in run OverDense_ExtremeRad. The mass fraction 
of ionized gas between 1-3 kpc from the galactic center is en- 
hanced by a factor of -200, compared to the fiducial run, even 
40 Myr after the PISN. Within the same distance range, the 
IGM temperature is a factor of -10 higher at this epoch. The 
relatively hot and ionized IGM in turn could affect accretion 
rates onto galaxies at later times. 

The effect of making PlSNe ten times more powerful in 
the H Mpc box was more drastic, as this box sampled lower 
mass galaxies. The baryon fraction in the first galaxy dropped 
to less than 10"^ after the first PISN compared to 1.7% in the 
standard run OverDense_nHle4_fid. The radii enclosing 80% 
of the baryons and metals are twice as large as in the standard 
run, demonstrating that the added energy in this extreme run 
coupled with the ISM more efficiently. Even with the standard 
feedback prescription we would expect a strong blowout in 



a halo of this mass (2.7 x 1O*M0 at this epoch). However, 
the extreme PISN energy is able to completely destroy the 
high-density gas clouds needed for star formation and prevent 
re-accumulation of dense gas for -100 Myr, when a merger 
brings in a fresh resevoir. 

These tests indicate that given enough energy input, the host 
halos of Pop 111 stars can become completely devoid of gas 
for cosmologically-significant intervals of time, particularly 
when they are below the mass threshold -3 x 1O^M0. How- 
ever, for the feedback parameters currently considered realis- 
tic (our fiducial runs), the feedback of the first stars is limited 
as illustrated in Figures |6] and |7] 

In the run with low masses of Pop 111 stars (Over- 
Dense_LowMass), without any PISN, metal transport is ex- 
tremely ineffective. At a = 0.055, 80% of the metals that have 
been generated by stars in the first star-forming galaxy are 
confined within 75 pc of the galactic center, compared to 550 
parsecs in the fiducial run. This test demonstrates that if Pop 
III stars did not have a top-heavy IMF, their contribution to en- 
riching the IGM with metals would be further marginalized. 

4. DISCUSSION AND CONCLUSIONS 

We have presented the results of simulations that im- 
plemented primordial star formation in the cosmological 
code ART. We find that the effects of stellar feedback on 
the amount of baryons and metals within the first galax- 
ies depends strongly on galaxy mass. For the lowest- 
mass galaxies (M/, -10^ M^) our re sults are similar to those 
of Brommetal. (2003); Whalen et al. (2008); Wise & Abell 
(2008); Wise et al. (2012), with gas and metals often being 
driven well beyond the virial radius of the Pop III star's host 
galaxy. For more massive galaxies (M/, > 10^ M0), however, 
a single PISN is not effective in evacuating the galactic ISM, 
as suggested in Wise & Cen (2009). Feedback from Pop III 
stars does not typically inject enough energy into the massive 
halos to permanently photo-evaporate the gas, drive metal- 
rich outflows past the virial radius, or have any other long- 
term impact. While Pop III stars can temporarily expel gas 
and quench star formation, the ISM begins to replenish soon 
after the SN explosion, as accretion from filaments at this 
epoch is very fast. All galaxies considered in our analysis ap- 
pear to have more than 50% of the universal baryon fraction 
restored 100 Myr after the first Pop III supernova event. Met- 
als are ejected anisotropically, and can travel relatively longer 
distances through the diffuse IGM in directions normal to the 
dense filaments which feed galactic accretion. This means 
that it typically takes more time for the ejected metals to be 
re-accreted into the galaxy, but we have demonstrated that this 
re-accretion does indeed occur, even in low-mass galaxies. 

A dividing line of M/, ?a 3 x 10^ Mq is of apparent impor- 
tance in determining whether the energy injection from the 
supernova at the end of the star's life can expel gas and metals 
out to a significant distance. The concept of a dividing line 
between early galaxies which suffer from significant blowout 
from those which do not has been consid ered in prior work 
(e.g. ICiardi et aL 2000; Ri cotti et al1l2002h . However, our re- 
sults point to a considerably lower threshold than what had 
been expected, as all but the very first galaxies are apparently 
robust to PISN feedback. The strength of this conclusion is 
bolstered by our use of a very strong feedback model for Pop 
III stars (even in our fiducial runs). Rea listically, the firs t 
stars may have a lower characteristic mass (iGreifet al.ll2011h . 
which in turn would mean that PlSNe are less frequent. 
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In order for simulations to capture the full range of relevant 
effects from Pop III star formation, resolving halos around 
IO^Mq with a sufficiently large number of particles is criti- 
cal. With insufficient resolution (less than 1000 DM particles 
for 10'' M© halos), all galaxies seem to reach M > 10^ 
without having yet formed a star. Since these galaxies are al- 
ready beyond the M/j « 3 x 10'' Mq dividing line, few disrup- 
tive effects from Pop III feedback are observed in these sim- 
ulations. Aggressive super-Lagrangian refinement may help 
resolve star formation in halos of lower mass, but requires a 
prohibitively large number of computations. A more practical 
approach is to begin simulations with sufficiently high resolu- 
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